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Introduction 

This  research  project  concerns  the  development  of  MRI  arterial  spin  tagging  to  non-invasively 
measure  breast  tissue  perfusion.  MR  dynamic  first-pass  contrast-enhanced  imaging  has  shown 
that  malignant  and  benign  breast  lesions  can  be  distinguished.  However,  it  may  have  limited 
importance  in  clinical  breast  diagnosis  due  to  significant  false-negatives  and  false-positives.  The 
arterial  spin  tagging  technique  was  developed  to  measure  tissue  perfusion  parameters  without  the 
use  of  contrast,  and  has  been  successfully  demonstrated  in  brain  and  kidney.  The  specific  aims 
are  to  (1)  refine  arterial  spin  tagging  pulse  sequences  and  imaging  protocols,  (2)  develop 
automated  data  analysis  software  for  measurement  of  breast  tissue  parameters,  and  (3)  compare 
the  technique  to  first-pass  contrast-enhanced  MRI  and  biopsy.  We  will  test  the  hypothesis  that 
arterial  spin  tagging  provides  accurate  and  precise  discrimination  between  normal  tissue,  benign 
and  malignant  lesions,  when  differences  in  perfusion  and  7j  exist.  Lesions  will  have  been 
previously  detected  by  clinically  accepted  diagnostic  imaging  procedures,  and  by  biopsy. 
Statistical  analysis  will  be  performed  to  access  the  correspondence  between  arterial  spin-tagging 
and  biopsy,  and  to  establish  the  relative  value  of  spin  tagging  compared  to  first-pass  contrast- 
enhanced  MRI.  We  hope  to  establish  that,  relative  to  using  first-pass  contrast-enhanced  imaging, 
false  positives  and  negatives  are  reduced  using  arterial  spin  tagging  by  virtue  of  increased  image 
signal-to-noise  ratio  (SNR),  higher  spatial  resolution,  and  the  unique  ability  to  obtain  estimates 
of  macromolecular  bound  fluid  fractions.  The  scope  of  effort  on  the  project  is  mainly  limited  to 
the  technical  aspects  of  development  of  a  new  methodology.  However,  the  project  also  includes  a 
rigorous  performance  comparison  with  the  current  gold-standard  methodology.  The  technique 
will  be  evaluated  in  sixty  patients,  with  roughly  equal  numbers  of  benign  and  malignant  lesions. 
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Body 

Description  and  Chronology  of  Project  Schedule  Changes 
Summary  of  Project  delays 

The  major  deficiency  in  our  progress  is  with  the  third  technical  objective.  Specifically, 
comparison  of  the  arterial  spin  tagging  sequence  with  the  dynamic  contrast  enhanced  first  pass 
study  has  not  yet  been  achieved.  Roughly  60  patients  should  have  been  scanned  at  this  point  in 
the  project,  but  only  two  have  been  scanned.  This  deficit  is  due  to  the  delay  in  finding  a 
replacement  physician  for  the  co-I  on  the  project,  Rebecca  Zulim,  who  left  UCDMC  November 
1,  1998.  Also,  UCD  Medical  Center  Radiology  Department  upgraded  all  their  MRI  systems  from 
Signa  Genesis  5.x  to  Signa  LX  8.2.5,  resulting  in  several  months  during  which  our  pulse 
sequences  were  not  compatible  with  any  of  the  MRI  systems  we  had  access  to.  Substantial  time 
was  needed  to  convert  the  pulse  sequences  to  the  LX  8.2.5  platform  from  the  Genesis  5.x 
platform.  This  time  period  for  upgrading  the  UCD  Med  Center  MRI  systems  was  not  planned  in 
advance,  and  the  time  to  convert  our  research  software  was  not  written  into  the  approved  Work 
Statement.  The  graduate  student  on  this  project,  David  Zhu,  and  I  have  worked  on  the  conversion 
of  the  arterial  spin  tagging  pulse  sequences,  and  on  the  development  of  new  arterial  spin  tagging 
sequences  based  on  our  “odd-hybrid”  EPI  technique  [Buonocore  MH,  Zhu  DC.  High  spatial 
resolution  EPI  using  an  odd  number  of  interleaves.  Magnetic  Resonance  in  Medicine  41  (6): 

1 199-1205  (1999)].  Also,  advances  have  been  made  in  our  image  processing  software  for  the 
project.  Because  of  the  delays,  the  budget  and  schedule  was  officially  pushed  back  6  months, 
into  the  3rd  year  which  had  a  light  work  schedule.  For  full  details,  please  see  the  referenced  letter 
dated  April  22,  1999,  and  the  approved  Budget  change,  in  the  appendix  material.  The  project 
now  continues  into  the  3rd  year,  as  dictated  by  approved  revised  schedule  and  budget.  Dr.  PD 
Schneider  is  the  co-investigator  on  the  project. 

Project  timeline  changes 

Rearrangement  of  budget  and  project  schedule 

April  1999:  With  the  realization  that  original  co-investigator  Rebecca  Zulim,  MD,  having 
moved  to  Kaiser  Permanente,  was  no  longer  going  to  be  effective  as  co-investigator  on  this 
project,  I  submitted,  on  April  22,  1999,  a  long  letter  explaining  the  situation,  along  with  a 
revised  budget  to  the  US  Army-BCRP  office.  Because  the  work  schedule  in  the  third  year 
was  relatively  light,  I  requested  transfer  of  remaining  salary  support  and  unexpended  funds 
from  the  2nd  year,  to  the  third  year  starting  Oct  1, 1999.  Over  half  of  the  salary  support 
was  unexpended  (March  1999  through  October  1999  salary  was  transferred).  Supplies 
costs  were  also  carried  forward.  The  stipend  and  tuition  remission  of  the  graduate  assistant, 
David  Zhu,  could  not  be  deferred.  He  continued  to  work  toward  the  specific  aims  of  the 
project  throughout  the  project  period.  The  letter  dated  April  22,  1999,  and  the  approved 
Budget  change,  are  included  in  the  appendix  material  of  this  Report. 

May  1999:  The  revised  budget  and  schedule  is  approved  by  the  US  Army  BCRP. 

Major  events  in  Co-Investigator  change 

Nov  1, 1998:  Rebecca  Zulim  quits  UCD  Medical  Center,  and  begins  work  at  Kaiser 
Permanente  in  Sacramento.  There  is  no  replacement  for  her  at  UCD  Medical  Center. 

Rebecca  did  not  inform  me  of  her  intent  to  leave  UCD.  Apparently,  the  Surgical  Oncology 
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Department  was  not  aware  of  it  either,  since  there  was  no  recruitment  effort  in  place  at  that 
time.  At  this  point,  I  had  no  physician  on  the  project  to  recruit  patients. 

December  1999:  Rebecca  expresses  interest  in  continuing  the  research  collaboration  while 
she  is  at  Kaiser.  As  PI,  I  rewrite  the  Protocols  and  Consent  Form  for  recruitment  of 
patients  at  Kaiser,  and  submit  these  to  our  Human  Subjects  Review  Committee  (HSRC). 
Rebecca  submits  same  protocol  and  consent  to  the  Kaiser  HSRC. 

March  1999:  There  is  still  no  response  from  the  Kaiser  HSRC.  Rebecca  continues  to  say 
that  she  wants  to  collaborate,  but  it  is  becoming  more  clear  that  it  will  not  work  out.  I 
received  the  UC  Davis  HSRC  committee  action  on  the  revised  Protocol  and  Consent  Form 
that  I  submitted  in  order  to  recruit  patients  from  Kaiser.  The  UC  Davis  HSRC  committee 
requires  minor  corrections  for  approval  of  the  protocol,  but  I  do  not  submit  these 
corrections,  due  to  the  non-responsiveness  of  Kaiser’s  HSRC. 

April  1999: 1  contact  Phil  Schneider,  MD-PhD,  Director  of  the  Surgical  Oncology  at  UCD 
Medical  Center,  to  let  him  know  that  Rebecca  is  not  going  to  collaborate  at  Kaiser.  I 
request  his  participation,  or  that  of  a  clinical  fellow.  He  agrees  to  collaborate,  and  I  rewrite 
the  Protocol  and  Consent  Form  a  second  time  to  revert  back  to  patient  recruitment  at  the 
UCDMC  Oncology  Clinic,  and  to  include  him  as  collaborator. 

May  24, 1999:  the  UC  Davis  HSCR  approves  the  modified  protocol  and  consent, 
reestablishing  our  intent  to  recruit  from  UC  Davis  Clinics.  The  only  modification  from  the 
original  approved  consent  is  the  inclusion  of  Dr.  Schneider  (and  removal  of  Dr.  Zulim).  I 
also  prepared  a  advertising  flyer  which  is  included  in  the  revised  consent. 

July  1999-Oct  1999:  Patients  are  not  being  recruited,  due,  according  to  Dr.  Schneider,  to 
his  preoccupation  with  other  projects.  However,  he  reconfirms  his  commitment  to  the 
project  and  suggests  that  I  visit  the  clinic  each  week  to  post  the  flyer. 

Major  events  in  MR1  system  upgrades 

December  1998:  The  main  hospital  MRI  is  upgraded  to  Horizon  LX  Echospeed,  Version 
8.2,  from  Signa  Advantage  OS  Version  5.4.  This  invalidates  our  pulse  sequence  software. 
Because  the  Imaging  Center  MRI  is  not  yet  available  for  research  use,  there  is  no  system  at 
UC  Davis  on  which  patient  studies  can  be  done. 

December  17, 1998:  Purchase  order  is  placed  with  GE  for  upgrade  of  the  MRI  Signa 
Advantage  System  at  the  UCD  Imaging  Center  to  a  high-performance  CV/i  system.  The 
intent  of  the  Radiology  Dept,  is  to  convert  this  system  into  one  dedicated  to  MRI  research 
The  Radiology  Department  Chair  designates  the  Imaging  Center  system  as  the  official 
MRI  Research  System.  All  Breast  Imaging  Studies  are  to  be  done  at  this  site,  as  soon  as  the 
system  is  delivered. 

December  1998-April  1999:  The  graduate  student  and  I  spend  our  time  on  the  project 
learning  the  EPIC  pulse  sequence  for  LX  environment,  obtaining  relevant  source  codes  and 
programs  for  conversion  of  research  pulse  sequences.  Acquisition  of  the  EPIC  8.2.5 
compilers  from  GE  was  delayed  due  to  our  research  agreement  negotiations,  which  were 
completed  in  May  1999.  We  currently  have  a  five-year  agreement  starting  May  1,  1999. 
Acquisition  of  the  EPIC  software  was  also  delayed  due  to  GE  administrative  delays.  For 
example,  the  EPIC  compilers  for  Version  8.2.5  were  not  obtained  until  issues  regarding  the 
long-term  research  agreement  were  resolved 

April  1999:  Research  MRI  is  installed  at  the  Imaging  Center.  Work  begins  on  conversion 
of  pulse  sequences.  Significant  performance  problems  with  the  system  delay  getting  the 
software  converted  and  running  properly  for  use  in  patients.  Full  payment  for  the  system  is 


not  made.  90%  of  payment  made  in  July  1999,  remaining  10%  to  be  made  in  November, 

1999,  with  stipulations. 

April  1999-Sept  1999:  The  arterial  spin  tagging  breast  sequence  based  is  converted  from 
Version  5.4  to  8.2.5,  and  uses  the  enhanced  gradients  of  the  CV/I  high  performance 
system.  In  addition,  the  odd-number  interleaf  EPI  sequence  (see  appendix  for  copy  of  this 
publication)  is  converted.  Current  work  involves  incorporating  this  acquisition  into  the 
arterial  spin  tagging  sequence.  In  other  words,  instead  of  SPGR  acquisition  with  short  TR, 
there  is  EPI  acquisition.  We  test  both  the  “odd-number  interleaf’  EPI  sequence,  and  the 
even-number  interleaf  sequences  that  GE  provided  through  the  UC  Davis  Research 
Agreement.  For  full  details,  see  section  below  on  pulse  sequence  development. 

Discussion  of  Approved  Statement  of  Work 

Technical  objective  1 

Task  1;  Months  1-6:  Implementation  and  testing  of  magnetization  transfer  pulses  for 
both  arterial  tagging  and  first-pass  contrast  enhanced  sequences. 

Report  for  1999 

Magnetization  transfer  (MT)  has  been  implemented  with  the  EPI-based  spin  tagging  sequence 
(see  Figure  1).  The  MT  radiofrequency  (RF)  pulse  selectively  excites  the  protons  of  the 
macromolecules  and  its  hydration  layer  by  a  90°  flip  angle.  Two  gradient  spoiling  pulses  (in  the 
x-  and  y-  directions)  are  applied  immediately  after  the  RF  pulse  to  spoil  the  transverse 
magnetization.  The  usual  RF  excitation  pulse  for  data  acquisition  is  applied  immediately  after 
these  spoiling  pulses.  Due  to  hardware  and  or  compiler  problems  that  are  currently  being 
investigated  by  GE,  this  pulse  sequence  has  been  implemented  only  on  pulse  sequence 
simulation  software  (EPIC),  not  yet  on  the  MRI  system.  The  hardware  problems  are  somehow 
preventing  this  sequence  from  running  successfully  on  the  MRI  system.  Testing  of  MT  using  the 
fast  SPGR-based  sequence  revealed  that  it  would  not  be  effective.  Because  the  magnetization 
transfer  RF  pulse  is  long  duration  (16  ms),  the  “small  TR”  (optr)  period  of  the  fast  SPGR 
sequence  was  effectively  doubled  in  length,  thus  sacrificing  short  scan  time  and  distorting  the 
signal  upon  which  the  perfusion  measurement  is  based. 

From  1998  Annual  Report  ... 

Magnetization  transfer  pulses  have  not  been  implemented  yet.  Mainly,  this 
delay  was  due  to  the  fact  that  the  MRI  system  at  UC  Davis  Medical  Center  is 
scheduled  to  be  upgraded  to  the  new  Horizon  LX  system  soon,  and  we  opted  to 
defer  this  pulse  sequence  development  until  we  were  using  this  new  LX 
software  platform.  There  have  been  multiple  delays,  but  this  upgrade  is 
scheduled  to  occur  by  January  1999.  The  existing  pulse  sequence  will  be 
converted  from  the  5.4  OS  platform  to  the  LX2  platform  (a  substantial 
change) ,  and  the  magnetization  transfer  pulses  will  be  implemented  also. 

Because  of  concerns  regarding  the  ultimate  utility  of  the  magnetization 
transfer  technique,  we  have  decided  that  implementation  and  testing  of 
arterial  spin- tagging  sequence  based  on  echo  planar  imaging  (EPI)  data 
acquisition  is  a  higher  priority.  EPI  based  spin-tagging  technique  can 
theoretically  improve  the  accuracy  of  the  Ti  estimation,  which  we  have  studied 
and  worked  with  extensively  (see  below) .  The  major  problem  with  magnetization 
transfer  is  that  it  represents  a  perturbation  on  an  already  small  signal. 
Therefore,  whether  we  will  observe  a  magnetization  transfer  effect  is 
questionable.  Nevertheless,  it  will  be  implemented  as  part  of  the  project. 
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Task  2:  Months  1-9:  Implementation  and  testing  of  interleaved  high-resolution  imaging 
technique,  for  both  arterial  tagging  and  first-pass  contrast  enhanced  sequences. 

Report  for  1999 

Summary 

Since  the  completion  of  the  upgrade  of  the  UC  Davis  Imaging  Center  MRI  system,  and  since  it’s 
designation  as  the  Research  MRI  system  for  UC  Davis,  we  have  thoroughly  studied  the  EPIC  LX 
8.2.5  programming  language  for  pulse  sequence  development,  and  the  operation  of  this  new  MRI 
system.  The  transition  to  the  new  research  MRI  required  a  complete  rewrite  of  all  pulse  sequence 
software  that  we  had  developed  for  this  project.  In  addition,  we  revised  the  pulse  sequences  to 
utilize  the  high  performance  gradients  (40  mT/m  peak,  150  mT/m/ms  rise)  of  the  system.  We 
believe  these  revised  sequences  will  provide  much  better  arterial  spin  tagged  data  for  perfusion 
measurement.  All  of  the  other  sequences  in  our  breast  imaging  protocol  were  set  up  for  the  new 
LX  8.2.5  platform  as  well. 

Conversion  of  the  Fast  SPGR  Based  Arterial  Spin  Tagging  Sequence 
The  authors  converted  the  fast  SPGR-based  arterial  spin  tagging  pulse  sequence  from  the 
Genesis  5.4  platform  to  the  LX  8.2.5  platform.  The  converted  sequence  was  also  improved  to 
achieve  higher  resolution  with  a  shorter  overall  scan  time.  The  new  MR  system  has  higher  peak 
gradient  amplitude  and  slew  rate.  Fast  data  acquisition,  which  requires  the  higher  read-out 
gradient  amplitude  (4.0  g/cm)  could  be  specified,  and  consequently  the  data  acquisition  period 
(optr)  was  shorten  from  16.62  ms  to  12.5  ms.  Since  optr  was  shortened,  a  larger  matrix  in  the 
phase  encode  direction  could  be  specified.  The  resolution  is  now  256  x  256  matrix  using  a  one- 
half  phase  field  of  view  acquisition.  Previously  we  used  a  256  x  240  matrix  with  one-half  phase 
field  of  view  acquisition.  The  total  time  for  the  image  acquisition  (the  TR  period),  including  the 
global  inversion  RF  pulse,  has  been  decreased  from  2.7  seconds  to  2.22  seconds.  The 
performance  of  this  new  sequence  has  been  confirmed  using  phantoms  and  one  test  subject. 

Development  of  an  EPI-Based  Arterial  Spin  Tagging  Sequence 
The  EPI-based  arterial  spin  tagging  sequence  was  rewritten  and  verified  on  the  research  MRI 
system.  The  data  acquisition  scheme  is  based  on  our  recently  published  high-resolution  odd- 
number  interleaf  EPI  sequence  [Buonocore  MH,  Zhu  DC.  High  spatial  resolution  EPI  using  an 
odd  number  of  interleaves.  Magnetic  Resonance  in  Medicine  41  (6):  1199-1205  (1999)].  This 
sequence  was  first  developed  and  verified  on  the  Signa  Advantage  1 .5  T  GE  MR  system. 
Unfortunately,  under  the  new  LX  system,  the  polarity  alternation  scheme  for  each  successive 
interleaf  we  had  implemented  caused  the  system  to  hang,  and  substantial  time  was  expended 
understanding  the  problem.  In  October  1999,  we  formulated  an  alternative  approach  to  eliminate 
the  requirement  of  polarity  alternation,  which  was  successfully  implemented.  One  persistent 
problem  that  we  have  with  all  interleaved  EPI  sequences  on  the  LX  system,  is  reduced  SNR 
relative  to  theoretical  predictions.  The  cause  of  this  is  not  known. 

In  the  new  arterial  spin  tagging  sequence  using  the  interleaved  EPI  acquisition,  the  following 
steps  are  performed  at  each  interleaf:  (1)  non-selective  RF  inversion  pulse,  (2)  20  EPI 
acquisitions,  each  of  which  is  preceded  by  a  10°  RF  pulse  (Figure  2).  These  20  acquisitions  are 
played  out  with  a  repetition  time  ( optr  in  Figure  2)  of  100  ms.  Each  acquisition  defines  one  point 
on  a  T\  recovery  curve.  The  total  time  for  each  interleaf  (TR  in  Figure  5.10)  is  2  seconds.  Data 
for  the  other  interleaves  is  collected  in  the  same  fashion  in  subsequent  TR  intervals.  Figure  3 
shows  the  loop  nesting.  The  completed  set  of  images,  showing  T1  recovery  of  the  magnetization, 
is  acquired  in  two  seconds  times  the  number  of  interleaves. 
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From  1998  Annual  Report ... 

An  interleaved  arterial  spin- tagging  technique,  that  allows  high  resolution 
(e.g.  256  x  256  or  more)  has  been  implemented  and  tested.  We  also  implemented 
a  rectangular  data  acquisition  technique,  which  acquires  data  sets  that  have 
different  numbers  of  points  along  the  frequency  and  phase  encode  direction. 
The  256  x  240  interleaved  acquisition  provides  the  best  spatial  versus 
temporal  resolution  tradeoff,  and  has  been  used  for  all  of  our  recent 
studies . 

Task  3:  Months  1-9:  Analysis  of  spin  tagging  sequence  to  understand  causes  of  existing 
baseline  offsets,  effects  of  inversion  slice  transition  profiles,  and  effects  of  RF  flip  angle 
profile  on  the  measurements,  with  implementation  and  testing  pulse  sequence 
modifications  to  minimize  these  imperfections. 

Report  for  1999 

Regarding  the  spin  tagging  sequence  utilizing  fast  SPGR,  these  issues  were  treated  in  the  1998 
report.  The  EPI  arterial  spin  tagging  sequence  has  not  been  evaluated,  but  the  same  solutions  for 
the  inversion  slice  transition  profiles,  and  effects  of  RF  flip  angle  profile,  will  apply. 

From  1998  Annual  Report  ... 

Matching  of  inversion  and  excitation  slice  profiles  was  greatly  improved  by 
using  customized  RF  pulses  designed  with  the  Shinnar-Le  Roux  (SLR)  algorithm. 
The  new  data  processing  technique  (based  on  a  semi-log  linear  regression  of 
inversion  time  (TI)  dependent  signal)  significantly  reduced  the  previously 
reported  problem  in  measuring  the  longitudinal  relaxation  time  (Tl)  caused  by 
baseline  offsets,  and  by  the  uncertainty  in  the  effective  inversion  time.  The 
optimal  RF  flip  angle  for  spin  excitation  during  spatial  encoding  was  found 
to  be  10  degrees,  based  on  experiments  over  a  range  of  flip  angles. 


Technical  objective  2 

Task  4:  Months  3-15:  Software  for  automatic  estimation  and  error  analysis  of 
perfusion,  tissue  water  longitudinal  relaxation  time,  and  extracellular  fluid  volume 
fraction  from  mathematical  models  and  user-defined  ROIs  from  spin  tagging  timeseries. 
Implementation  of  pharmacokinetic  model  calculations  and  error  analysis  for  first-pass 
contrast  enhanced  imaging. 

Report  for  1999 

From  October  1998  through  April  1999  (during  which  time  there  was  no  MRI  system  for  us  to 
develop  on,  or  do  patient  studies),  image  processing  software  development  continued. 
Improvements  in  the  BreastView  program,  development  of  “dispAlls”  program  for  analysis  of 
general  MRI  images,  and  implementation  of  statistical  analysis  of  breast  lesions  based  on 
Bayesian  approach,  was  completed.  Statistical  analysis  of  available  breast  lesion  data  continued 
(See  below).  The  most  important  new  function  of  BreastView  is  that  for  analysis  of  the  first-pass 
contrast  enhanced  images.  The  program  creates  an  image  based  upon  the  time  derivative  of  the 
rate  of  rise  of  the  signal,  so  that  the  contrast  enhancement  effect  can  be  better  visualized  (See 
Figures  4  and  5). 

From  1998  Annual  Report  ... 


A  software  program  (BrView)  has  been  written  in  C,  X  Window  System,  and 
Motif,  and  implemented  on  an  SGI  02  computer  (paid  for  by  grant  funds)  for 
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analyzing  and  visualizing  the  MR  images  and  timeseries.  Much  effort  has  been 
expended  on  developing  the  capability  to  easily  review  all  of  the  images 
taken  in  each  study,  and  cross-reference  pixel  locations.  Much  more  effort 
than  anticipated  has  been  expended  in  the  development  of  a  robust  technique 
for  identifying  so-called  "suspicious  regions"  based  on  the  Ti  and  perfusion 
measurements.  The  measurement  of  Ti,  upon  which  the  perfusion  measurement 
entirely  depends,  is  based  on  a  semi-log  linear  regression  technique 
developed  by  the  investigators.  The  Ti,  perfusion  (f/A)  and  standard  errors  of 
these  quantities  are  estimated  automatically  using  this  robust  technique.  The 
so-called  feature  images  display  these  quantities  in  color  and  grayscale. 

These  quantities  are  used  to  calculate  a  "suspicion  index"  for  each  pixel, 
and  thereby  identify  regions  of  breast  tissue  suspicious  for  malignancy. 

The  implementation  of  the  pharmacokinetic  model  calculations  and  error 
analysis  for  first-pass  contrast  enhanced  imaging  has  not  been  completed. 

Thus  far,  the  software  allows  the  user  to  click  on  the  reference  image  of  the 
dynamic  study,  to  display  the  time  profile  of  the  signal  at  that  pixel.  The 
implementation  of  the  pharmacokinetic  model  is  currently  a  major  focus  of  the 
current  effort. 

Task  5:  Months  3-15:  Software  for  registration  (including  implementation  and  testing 
of  motion  correction  and  physiological  noise  reduction  algorithms)  and  overlay  of 
images  from  high  resolution  T ;  weighted,  spin-tagging,  and  contrast  enhanced  studies. 

Report  for  1999 

First  described  in  the  1998  Annual  Report,  the  BreastView  program  has  been  continually 
upgraded  and  enhanced  with  the  latest  algorithms  for  processing  and  analysis  of  breast  imaging 
studies,  including  those  for  spin-tagged  and  dynamic  contrast  enhanced  images.  The  specific 
developments  this  year  of  BreastView,  and  of  another  program,  dispAlls,  are  described. 

BreastView 

Between  October  1998  and  April  1999,  the  BreastView  program  was  made  “clinician-friendly”. 
The  graphical  user  interface  now  guides  the  user  through  the  entire  visualization  and  analysis 
process  for  a  breast  imaging  study.  A  "Help"  button,  linked  to  extensive  set  of  manual  pages, 
now  appears  on  the  first  page  of  the  interface.  Clicking  the  Help  button  brings  up  an  html  file 
with  necessary  operating  instructions.  The  help  files  were  written  for  clinicians.  The  clinician 
can  now  begin  using  the  program  without  first  studying  the  written  instructions.  If  an 
inappropriate  option  is  selected,  the  program  now  provides  help  to  the  user  to  select  a  more 
appropriate  option  or  operation.  Several  new  features  and  functions  have  been  implemented,  as 
indicated  on  Figures  6  and  7.  The  most  important  new  function  of  BreastView  is  the  analysis  of 
the  first-pass  contrast  enhanced  images.  The  program  creates  a  contrast  enhancement- weighted 
image  based  upon  the  time  derivative  of  the  rate  of  rise  of  the  signal,  so  that  the  contrast  effect 
can  be  better  visualized  (See  Figures  4  and  5). 

The  dispAlls  program 

This  program  has  been  written  for  quickly  viewing  and  performing  ROI  analysis  on  one  or  a 
series  of  images.  The  program  is  particularly  useful  for  quickly  reviewing  the  timeseries  of 
signal  changes  in  an  ROI  from  the  dynamic  contrast  enhanced  studies.  The  program  displays  one 
image,  or  a  sequence  of  images,  with  or  without  an  image  header.  It  performs  ROI  analysis  and 
plots  the  calculated  values. 
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From  1998  Annual  Report  . . . 

Using  the  BrView  software  program,  image  pixels  containing  high  values  of 
perfusion  (f/X),  and  moderate  to  high  values  of  Ti,  and  low  standard  errors, 
are  identified  as  suspicious,  and  assigned  a  "suspicion  index"  based  on 
finding  similar  abnormal  values  at  spatially  adjacent  pixels.  Suspicious 
pixels  can  be  overlaid  on  the  high-resolution  Ti  and  Ti  clinical  images  by 
color  mapping,  and  simultaneously  presented  with  the  first-pass  contrast 
enhanced  timeseries  at  these  pixels.  This  provides  the  user  with  a 
comprehensive  anatomical  and  functional  view  of  the  suspicious  regions,  and 
facilitates  making  a  decision  regarding  the  malignant  nature  of  the  lesion. 
Finally,  a  motion  artifact  estimation  and  correction  algorithm  has  been 
developed  and  implemented  in  a  separate  software  program. 


Technical  objective  3 

Task  6:  Months  9-24:  60  patients  with  malignant  and  benign  breast  lesions  will  be 
imaged  using  T/  weighted,  first-pass  contrast  enhanced,  and  arterial  spin  tagging  MRI 
pulse  sequences. 

Report  for  1999 

Unfortunately,  no  additional  subjects  have  been  done.  This  is  due  to  the  loss  of  Rebecca  Zulim, 
MD,  as  the  co-investigator,  and  the  transition  to  a  dedicated  research  MRI  system  that  operated 
at  LX  8.2.5,  as  described  above.  Through  Sept  1998,  eighteen  subjects,  including  1 1  patients 
referred  by  their  physicians  and  seven  volunteers  referred  by  themselves,  had  participated  in  the 
breast  studies.  Three  patients  and  two  volunteers  were  eliminated  from  the  analysis  because  the 
studies  were  not  done  according  to  protocol,  due  to  excess  motion,  or  the  abnormal  mass  being 
too  close  to  the  chest  wall  to  be  measured  reliably,  or  no  patient  biopsy  to  confirm  the  MR  study 
result.  There  are  total  of  13  useful  cases  left  for  analysis.  For  all  these  cases,  the  suspicious  pixels 
were  identified  based  on  the  following  criteria:  7jn:  >  0.5  sec  f/X:  >  0.1  sec'1,  STD  of f/X:  <0.1 
sec'1,  and  the  suspicion  level  threshold  was  set  at  20.2%.  Analysis  has  been  done  using  the 
program  BreastView.  Table  8.5  shows  the  result  of  these  13  cases. 


Table  8.5.  Summary  of  Breast  Case  Studies 


Date  of  Subject 
Study  Referred 
by 

Palpable  Seen  on 
(Yes  or  Clinical 
No)  ?  MRI  ? 

Level  of 
Suspicion 
Based  on 
Tin  and  f/X 

Assessment 
Based  on 
Suspicion 
Level 

Abnormal  Mass 
Type 

05/03/97  Physician 

No 

Yes 

Low 

Negative 

Benign 

(Fibroadenoma) 

05/23/97  Physician 

Yes 

Indeter¬ 

minable 

Moderate 

High 

Positive 

Malignant 

(Ductal 

Carcinoma) 

08/01/97  Physician 

Yes 

No 

High 

Positive 

Malignant 

(Ductal 

Carcinoma) 

11/15/97  Physician 

Yes 

No 

Low 

Negative 

Benign 

12/10/97  Physician 

Yes 

No 

High 

Positive 

Benign 
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04/13/98  Physician 

Yes 

Indeter¬ 

minate 

Moderate 

High 

Positive 

Benign 

07/27/98  Physician 

Yes 

No 

Low 

Negative 

Benign 

08/07/98  Physician 

Yes 

No 

High 

Positive 

Benign 

02/07/97  Self 

No 

No 

Low 

Negative 

Normal 

03/05/97  Self 

Yes 

Yes 

Low 

Negative 

Benign 

(Fibroadenoma) 

07/26/97  Self 

No 

No 

Low 

Negative 

Normal 

11/16/97  Self 

Yes 

Yes 

Low 

Negative 

Benign 

(Fibroadenoma) 

12/14/97  Self 

No 

No 

Low 

Negative 

Normal 

Normal  =  normal  subject  who  does  not  have  any  tumor  based  on  reliable  health  history, 
Abnormal  Mass  Type  =  determined  based  on  biopsy  or  reliable  health  history.  Based  on  the 
summary  shown  in  Table  8.5,  the  number  of  true  positive  (TP)  is  2;  the  number  of  false  negative 
(FN)  is  0;  the  number  of  false  positive  (FP)  is  3  and  the  number  of  true  negative  (TN)  is  8. 
Therefore,  the  true-positive  fraction  or  sensitivity  is  100%,  calculated  by  TP/(TP+FN).  The  false¬ 
positive  fraction  is  27.3%,  calculated  as  FP/(FP+TN).  The  specificity  is  72.7%,  which  is  equal  to 
one  minus  the  false-positive  fraction.  Because  the  number  of  cases  investigated  is  small,  these 
figures  do  not  provide  sufficient  statistical  confidence  to  evaluate  the  technique.  After  more 
breast  cases  have  been  completed,  a  ROC  curve  can  be  generated  using  different  suspicion  level 
thresholds  or  different  thresholds  of  T\n,  ffk  and  STD  of ffk.  If  a  large  pool  of  case  studies  is 
available,  a  more  sophisticated  Bayesian  statistical  model  can  also  be  generated. 

From  1998  Annual  Report ... 

We  have  recruited  only  two  subject  studies  for  both  first-pass  contrast 
enhanced  and  arterial  spin  tagging  MRI  pulse  sequences.  Now  that  the  software 
development  is  substantially  completed,  we  are  trying  to  improve  the  rate  at 
which  patients  are  recruited  into  the  study.  It  is  likely  that  the  patient 
studies  will  extend  into  the  third  year  of  the  proposal.  We  are  currently 
discussing  with  physicians  at  a  local  Kaiser  Hospital  for  recruitment  of 
additional  patients.  Appropriate  documentation  will  be  submitted  for 
approvals  upon  reaching  a  firm  agreement. 

Task  7  (listed  as  Task  6  in  the  original  grant  application):  Months  12-36:  Automated 
pharmacokinetic  analysis,  blinded  image  reading,  and  statistical  comparison  of  arterial 
spin  tagging,  contrast  enhanced  MRI,  and  biopsy  results. 

Report  for  1999 

We  developed  an  analysis  program  for  the  detection  of  the  features  of  the  dynamic  contrast 
enhanced  study,  and  have  developed  a  statistical  approach  to  the  determination  of  the  malignant 
status  of  the  lesion. 

Bayesian  Statistical  Analysis 

The  identification  of  a  suspicious  pixel  based  upon  the  index  threshold,  that  was  developed  last 
year,  is  a  simple  approach  to  cancer  detection.  When  a  large  number  of  studies  becomes 
available,  a  Bayesian  decision  technique  [53]  will  be  applied.  This  technique  will  be  used  to 
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provide  ROC  (receiver  operating  characteristic)  curves  that  are  commonly  used  for  evaluation  of 
diagnostic  algorithms.  To  summarize  the  Bayesian  technique,  let  vector  space  x  =  ( T\,ffk ,  STD 
of f/X,  etc.)  for  a  ROI  (region  of  interest)  in  a  breast.  The  risks  associated  with  a  given  vector  x 
are 


Risk(Label  Benign  I  x)  =  r(Label  Benign  I  Malignant)  P(Malignant  I  x)  (7.1) 

+  r(Label  Benign  I  Benign)  P(Benign  I  x) 
and 

Risk(Label  Malignant  I  x)  =  r(Label  Malignant  I  Malignant)  P(Malignant  I  x)  (7.2) 

+  r(Label  Malignant  I  Benign)  P(Benign  I  x) 

where,  Risk(Label  Benign  I  x)  =  the  risk  value  to  label  a  ROI  "Benign"  for  a  given  x,  Risk(Label 
Malignant  I  x)  =  the  risk  value  to  label  a  ROI  "Malignant"  for  a  given  x,  r(Label  Benign  I 
Malignant)  =  the  risk  factor  to  label  a  ROI  "Benign"  but  actually  be  malignant,  r(Label  Benign  I 
Benign)  =  the  risk  factor  to  label  a  ROI  "Benign"  and  actually  be  benign,  r(Label  Malignant  I 
Malignant)  =  the  risk  factor  to  label  a  ROI  "Malignant"  and  actually  be  malignant,  r(Label 
Malignant  I  Benign)  =  the  risk  factor  to  label  a  ROI  "Malignant"  but  actually  be  benign, 
P(Malignant  I  x)  =  the  probability  of  being  malignant  for  a  given  x,  and  P(Benign  I  x)  =  the 
probability  of  being  benign  for  a  given  x.  Generally,  there  is  no  risk  for  correct  labeling,  and  thus 
r(Label  Benign  I  Benign)  =  0,  and  r(Label  Malignant  I  Malignant)  =  0.  Eq.  (7.1)  and  (7.2)  then 
become 

Risk(Label  Benign  I  x)  =  r(Label  Benign  I  Malignant)  P(Malignant  I  x)  (7.3) 

and 

Risk(Label  Malignant  I  x)  =  r(Label  Malignant  I  Benign)  P(Benign  I  x)  (7.4) 

Furthermore,  since, 

P(Malignantlx)  =  P<*'Malignant)P(Malignant)  (7.5) 

P(x) 


and 


P(Benignlx) 


p(xl  Benign)P(Benign) 

P(x) 


(7.6) 


we  have 


Risk  Ratio  =  Risk(Label  Benign  I  x) 

Risk(Label  Malignant  lx) 

_  r(Label  Benign  I  Malignant)  P(Malignant  I  x) 

r(Label  Malignant  I  Benign)  P(Benign  I  x)  (7.7) 

r(Label  Benign  I  Malignant)  p(x  I  Malignant)  P(Malignant) 
r(Label  Malignant  I  Benign)  p(x  I  Benign)  P(Benign) 
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.  r(Label  Benign  I  Malignant)  .  ..  .  ....  ^  p(x  I  Malignant) 

The  ratio  — - - - - -  is  estimated  by  the  clinician.  The  ratio  — - — — — ; — — 

r (Label  Malignant  I  Benign)  p(x  I  Benign) 

is  determined  from  a  training  data  set.  The  ratio  P^11^11^  can  be  obtained  from  literature  or 

P(Bemgn) 

estimated  from  the  training  data  set;  but  the  source  of  the  ratio  obtained  must  correspond  with  the 
source  of  subjects  used  in  the  breast  study.  If  the  Risk  Ratio  is  greater  than  1,  then  the  ROI 
would  be  labeled  "Malignant";  otherwise,  the  ROI  would  be  labeled  "Benign".  Or 


Risk  Ratio 


Label  Malignant 

S 

Label  Benign 


The  values  of  risk  factors  are  determined  by  the  clinician  or  the  researcher  based  on  his/her 
experience.  For  example,  if  the  clinician  thinks  that  it  is  too  risky  to  label  a  malignant  mass  to  be 
"Benign",  he/she  would  give  a  high  value  for  r(Label  Benign  I  Malignant),  and  a  relative  low 
value  for  r(Label  Malignant  I  Benign).  Choosing  the  different  values  of  the  risk  factor  ratio 

r(Label  Benign  I  Malignant)  eqU  jvajent  t0  choosing  different  values  of  the  decision  threshold. 

r(Label  Malignant  I  Benign) 

The  pairs  of  false  positive  and  true  positive  fractions  resulting  from  the  different  values  of  the 
risk  factor  ratio  would  allow  the  plotting  of  the  ROC  curve,  which  can  be  used  to  evaluate  the 
performance  of  the  technique  presented  [1]. 

™  .  p(x  I  Malignant)  ,  P(Malignant)  .  .  .  ,  .  .  .  , 

The  two  ratios,  — - - - and - - - ,  requires  a  large  training  data  set  to  be 

p(x  I  Benign)  P(Benign) 

statistically  significant.  This  requires  more  breast  case  studies  in  the  future. 
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MT  pulse 


Figure  1.  The  EPI  based  spin  tagging  sequence  with  MT  pulses  incorporated. 
The  odd-number  interleave  EPI  technique  discussed  is  implemented  in  data 
acquisition.  The  upper  figure  shows  the  timing  diagram  for  each  interleave  at 
acquisition.  A  180°  degree  RF  pulse  is  followed  by  N  repeated  EPI  data 
acquisitions  (as  shown  at  the  lower  figure). 


z 

gradient 


Figure  2.  The  EPI  based  spin  tagging  sequence.  The  odd-number  interleave  EPI 
technique  discussed  is  implemented  in  data  acquisition.  The  upper  figure  shows 
the  timing  diagram  for  each  interleave  acquisition.  An  180°  degree  RF  pulse  is 
followed  by  N  repeated  EPI  data  acquisitions  (as  shown  in  the  lower  figure). 
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Figure  3.  The  loop  nesting  for  the  EPI  based  spin  tagging  sequence. 
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Figure  6.  The  opening  page  of  the  BreastView  program. 


Figure  7.  Aug  7, 1998  breast  study.  Suspicion  level  map  is 
created  using  the  “Create  Detection  Image”  button  after 
setting  the  thresholds. 


Key  Research  Accomplishments 

A.  Pulse  sequences  that  are  compatible  with  a  Signa  Horizon  LX  EchoSpeed  CV/i  system,  running  LX 
8.2.5  system  software: 

1 .  Fast  SPGR-based  arterial  spin  tagging  sequence 

2.  Odd-number  interleaved  EPI  sequence 

3.  Interleaved  EPI-based  arterial  spin  tagging  sequence 

B.  Imaging  processing  software 

1.  BreastView 

2.  DispAlls 

3.  Bayesian  analysis  of  arterial  spin  tagged  image  data 

Reportable  Outcomes 

A.  Pulse  sequences  for  Version  8.2.5  (MH  Buonocore  and  D  Zhu): 

1.  Arterial  spin  tagging  sequence  with  Fast  SPGR  Acquisition 

2.  Odd-number  interleaf  EPI 

3.  Arterial  spin  tagging  sequence  with  Odd-number  interleaf  EPI  acquisition 

4.  Arterial  Spin  Tagging  sequence  with  Even-number  hybrid  EPI  acquisition 

B.  Breast  Imaging  Protocols  (MH  Buonocore  and  D  Zhu) 

1 .  Anatomical,  spin  tagged,  and  dynamic  c-e  breast  imaging  protocol 

C.  SGI  Graphics  Software  for  Breast  image  analysis  (MH  Buonocore  and  D  Zhu): 

1 .  BreastView  (display  program  for  breast  imaging  protocol) 

2.  DispAll  (display  program  for  general  MRI  images) 

D.  PhD  Degree  in  Biomedical  Engineering,  December  1999  (D  Zhu) 

E.  Employment  as  GE  Applications  Engineer,  starting  Jan  2000  (D  Zhu): 

Conclusions 

We  have  successfully  converted  and  enhanced  the  arterial  spin  tagging  pulse  sequences  to  work  on  a  new 
research  MRI  system,  and  have  made  several  enhancements.  Understanding  of  the  8.2.5  pulse  sequence 
language  on  the  new  system  represents  a  major  advance  of  the  project.  These  pulse  sequences  will  be  used 
in  all  future  work,  and  hopefully  there  will  not  be  a  substantial  MRI  system  upgrade  such  as  we 
experienced  for  a  few  years  at  least.  Image  display  and  statistical  analysis  software  has  been  written  for 
the  spin  tagging  and  contrast  enhanced  studies,  and  to  allow  analysis  to  be  done  by  a  busy  clinician 
without  computer  expertise. 
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High  Spatial  Resolution  EPI  Using 
an  Odd  Number  of  Interleaves 

Michael  H.  Buonocore*  and  David  C.  Zhu 


Ghost  artifacts  in  echoplanar  imaging  (EPI)  arise  from  phase 
errors  caused  by  differences  in  eddy  currents  and  gradient 
ramping  during  left-to-right  traversal  of  kx (forward  echo)  versus 
right-to-left  traversal  of  kx  (reverse  echo).  Reference  scans  do 
not  always  reduce  the  artifact  and  may  make  image  quality 
worse.  To  eliminate  the  need  for  reference  scans,  a  ghost 
artifact  reduction  technique  based  on  image  phase  correction 
was  developed,  in  which  phase  errors  are  directly  estimated 
from  images  reconstructed  separately  using  only  the  forward  or 
only  the  reverse  echos.  In  practice,  this  technique  is  applicable 
only  to  single-shot  EPI  that  produces  only  one  ghost  (shifted  % 
the  field  of  view  from  the  parent  image),  because  the  technique 
requires  that  the  ghosts  do  not  completely  overlap  the  parent 
image.  For  higher  spatial  resolution,  typically  an  even  number  of 
separate  k-space  traversals  (interleaves)  are  combined  to  pro¬ 
duce  one  large  data  set.  In  this  paper,  we  show  that  data 
obtained  from  an  even  number  of  interleaves  cannot  be  com¬ 
bined  to  produce  only  one  ghost,  and  image  phase  correction 
cannot  be  applied.  We  then  show  that  data  obtained  from  an  odd 
number  of  interleaves  can  be  combined  to  produce  only  one 
ghost,  and  image  phase  correction  can  be  applied  to  reduce 
ghost  intensity  significantly.  This  “odd-number  interleaf  EPI” 
provides  spatial  and  temporal  resolution  tradeoffs  that  are 
complementary  to,  or  can  replace,  those  of  even-number  inter¬ 
leaf  EPI.  Odd-number  interleaf  EPI  may  be  particularly  useful  for 
MR  systems  in  which  reference  scans  have  been  unreliable. 
Magn  Reson  Med  41:1199-1205, 1999.  ©  1999  Wiley-Liss,  Inc. 
Key  words:  magnetic  resonance  imaging;  interleaved  echo  pla¬ 
nar  imaging;  ghost  artifact  suppression;  image  phase  correc¬ 
tion 

To  increase  spatial  resolution  in  echoplanar  imaging  (EPI), 
two  or  more  separate  k-space  traversals  (interleaves)  can  be 
used  to  acquire  different  sets  of  k-space  lines  that  are 
combined  to  produce  one  large  data  set.  Typically,  this 
“interleaved  EPI”  uses  an  even  number  of  interleaves 
(1-4).  Interleaved  EPI  is  particularly  sensitive  to  the  genera¬ 
tion  of  ghost  artifacts,  which  arise  from  the  differences  in 
eddy  currents  and  gradient  ramping  that  occur  during 
left-to-right  traversal  of  kx  (forward  echo,  denoted  J)  versus 
right-to-left  traversal  of  kx  (reverse  echo,  denoted  /).  These 
forward  and  reverse  echos  cause  different  phase  errors  that 
result,  upon  two-dimensional  (2D)  Fourier  transformation, 
in  ghost  artifacts  in  the  phase  encode  (/cy)  direction.  These 
ghosts  can  often  be  reduced  using  phase  error  information 
obtained  from  non-ky  encoded  reference  scans  (5-11). 

Previously,  one  of  the  authors  published  a  ghost  artifact 
reduction  technique  (12)  based  on  extracting  the  spatially 


dependent  phase  errors  directly  from  images  reconstructed 
separately  using  only  the  forward  or  only  the  reverse 
echos.  This  technique  was  developed  only  for  single-shot 
EPI,  in  which  forward  and  reverse  echos  are  alternated  in 
k-space,  and  Fourier  reconstruction  produces  a  single 
ghost  artifact  shifted  by  V2  the  field  of  view  (FOV)  from  the 
location  of  the  “parent”  image.  The  technique  requires,  for 
each  location  in  the  frequency  encode  direction,  at  least 
one  pixel  within  the  parent  image  that  the  ghost  artifact 
does  not  overlap.  The  technique  could  not  be  applied  to 
even-number  interleaf  EPI  that  have  been  presented  in  the 
literature,  because  the  multiple  ghost  artifacts  produced  in 
these  scans  entirely  overlapped  the  parent  image.  With  h 
interleaves,  ghost  and  parent  images  were  separated  by 
only  ±l/2h  the  FOV. 

This  paper  introduces  and  illustrates  the  advantages  of 
using  an  odd  number  of  interleaves.  The  fundamental 
criteria  for  the  interleaved  EPI  sequence  design  was  to 
achieve  a  k-space  trajectory  that  mimics  the  following 
features  of  a  single-shot  EPI  method  that  are  critical  to 
producing  images  without  ghost  artifacts:  1)  smooth  T2* 
decay  of  the  signal  as  a  function  of  ky\  and  2)  phase  errors 
due  to  forward  and  reverse  echos  arranged  in  an  alternating 
/-r-/-r-7-r-/-r-etc.  pattern  as  a  function  of  ky.  Reduction  of 
ghost  artifacts  from  TJ  decay  has  been  previously  ad¬ 
dressed  in  the  literature  (2).  In  single-shot  EPI,  the  decay  is 
smooth  as  a  function  of  ky  and  causes  only  a  local  blurring 
of  the  image.  In  interleaved  EPI,  the  decay  can  significantly 
affect  image  quality  if  the  combined  data  set  is  strongly 
modulated  as  a  function  of  ky.  Strong  modulation  is 
avoided  by  temporally  delaying  each  interleaf  in  accor¬ 
dance  with  its  starting  ky,  such  that  the  combined  data  set 
acquires  a  smooth  T£  decay  in  the  ky  direction.  As  we  will 
show,  the  need  to  assemble  a  1-r-l-r-l-r-l-r-etc.  alternation  of 
phase  errors  in  the  interleaved  EPI  data  set  leads  naturally 
to  the  use  of  an  odd-number  of  interleaves.  An  odd  number 
is  necessary  to  generate  a  single  ghost  that  is  well  separated 
from  the  parent  image  and  to  utilize  image  phase  correc¬ 
tion. 


THEORY 

Let  N  equal  the  number  of  raw  data  points  in  the  kx  and  ky 
directions  and  in  the  final  image.  Define  gjr(n)  according  to 
the  direction  of  the  nth  k-space  line, 


_  1  if  nth  line  acquired  from 
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locations  and  complex  amplitude  factors  of  resulting  ghosts, 
define 


gi(ri)  =  lMl  +  gir(n))  and  gr(n)  =  '/2(1  -  g,r(n))  (2] 


as  functions  identifying,  respectively,  the  forward  and 
reverse  echos  as  a  function  of  k-space  line  indexed  by  n. 
Set  the  spatial  scale  such  that  the  object  being  recon¬ 
structed,  M(x,y),  is  defined  at  N  x  N  integer- valued  points 
(x,  y)  in  the  spatial  domain.  Following  from  Eq.  [3]  in  ref. 
12,  the  signal  at  points  labeled  by  integers  in  k-space  can 
then  be  written  as 


S(m,  n)  —  (1/A/2) 


N/2-1  N/  2-1 


2  2 

x^-N/2  y=-N/2 


M{x,  y) 


X  (gj(n)  exp  (i0  (x,  y)) 


[3] 


+  gr{n)  exp  (-J0(x,  y))) 


FIG.  1.  Three-interleaf  EPI  pulse  sequence  diagram.  For  each 
interleaf,  relative  temporal  delays  and  polarity  changes  of  x  (fre¬ 
quency  encode),  y  (phase  encode),  and  z  (slice-select)  gradient 
waveforms  are  shown  for  each  of  the  three  interleaves  (labeled  #1 , 
#2,  and  #3).  Interleaf  #1  acquires  the  1st  k-space  line  and  every  third 
thereafter.  Interleaves  #2  and  #3  acquire  the  2nd  and  3rd  k-space  lines 
and  every  third  line  thereafter,  but  the  first  acquired  k-space  line  from 
each  interleaf  is  discarded.  To  create  a  smooth  T2*  decay  in  the 
combined  data  set,  interleaves  #2  and  #3  are  temporally  shifted  such 
that  their  2nd  acquired  k-space  line  is  offset  by  1/3  and  2/3  of  the  total 
time  for  one  trapezoid  of  the  readout  gradient.  Phase  encoding 
occurs  at  the  transitions  between  forward  and  reverse  echos  and  is 
also  temporally  shifted  in  each  interleaf.  The  slice-select  gradient  has 
a  constant  spoiler  gradient  after  data  acquisition  to  reduce  residual 
transverse  magnetization.  Z-gradient  and  RF  pulses  for  slice  selec¬ 
tion  are  identical  for  all  interleaves. 


due  to  gradient-dependent  main  field  offsets,  echo  delays, 
and  echo  drifts.  If  the  data  alternate  between  forward  and 
reverse  echos  in  successive  k-space  lines,  then  gir(n)  =  1  for 
n  even,  and  gir[n)  —  - 1  for  n  odd.  This  alternating  pattern  of 
echos  produces  a  single  ghost  shifted  Vz  the  FOV,  which 
can  be  suppressed  using  image  phase  correction  (12). 

Multiple  ghost  artifacts,  some  shifted  only  a  small  dis¬ 
tance  from  the  parent  image,  occur  with  more  complicated 
patterns  of  forward  and  reverse  echos  as  a  function  of  ky.  To 
derive  the  relation  between  the  echo  pattern,  and  the 


Table  1 

Interleaved  EPI  Used  for  Figs.  4-7 


No.  of 
interleaves 

Data  matrix 
of  each  interleaf1 

Resulting 
image  matrix1 

1 

64  X  64 

64  x  64 

2 

128  x  64 

128  x  128 

3 

128  X  44 

128  x  128 

9 

256  X  30 

256  x  256 

Matrices  are  given  as  frequency  encode  x  phase  encode. 


2ir  /  2tt  \ 

*  exp  \-i  —  mx  exp  \~i  —  nyj 


Regrouping  terms,  and  carrying  out  an  Inverse  Fourier 
transform  on  Eq.  [3],  the  reconstructed  image,  M(x,  y),  can 


#3- 

#1- 


k: 

k; 


v 


‘  V  #2 


FIG.  2.  Odd-number  interleaf  EPI.  Three  separate  interleaves  (la¬ 
beled  #1,  #2,  and  #3)  are  combined  to  produce  a  l-r-l-r-l-r-l-r-etc. 
pattern  of  echos.  In  our  implementation  of  128  x  128  interleaved  EPI, 
each  interleaf  acquires  44  k-space  lines,  yielding  an  excess  of  4  lines 
that  are  discarded.  Interleaf  #1  acquires  the  1st,  4th,  7th,  etc.  k-space 
lines  (counting  from  1  at  top),  interleaf  #2  acquires  the  2nd,  5th,  8th, 
etc.  lines,  and  interleaf  #3  acquires  the  3rd,  6th,  9th,  etc.  lines. 
Interleaves  #1  and  #2  are  initiated  with  a  forward  echo,  while  interleaf 
#3  is  initiated  with  a  reverse  echo.  The  first  acquired  k-space  line  from 
interleaf  #2,  and  the  first  acquired  line  from  interleaf  #3,  are  not  used 
in  the  combined  data  set.  Also,  the  last  k-space  line  (131th  in 
combined  data  set)  from  interleaf  #3,  and  the  last  line  from  interleaf 
#1  (132th  in  combined  data  set)  are  discarded.  Discarding  the  first 
k-space  lines  is  based  on  the  presumption  that  equilibrium  phase 
errors  may  occur  only  after  the  first  line,  meaning  that  the  2nd  and  3rd 
lines  of  the  combined  data  set  should  not  be  acquired  as  the  1 st  line  of 
one  of  the  interleaves. 


Odd-Number  Interleaf  EPI 
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FIG.  3.  Even-number  interleaf  EPI.  In  our  implementation  of  128  x 
128  interleaved  EPI,  each  interleaf  acquires  32  k-space  lines,  and  no 
lines  are  discarded.  Four  separate  interleaves  (labeled  #1 ,  #2,  #3, 
and  #4)  are  combined  to  produce  different  echo  patterns.  Each  figure 
shows  exactly  two  cycles  of  the  smallest  repeating  pattern  of  the 
combined  data  set.  a:  A  l-l-l-l-r-r-r-r-etc.  pattern  is  formed  by  using 
interleaf  #j  (j  =  1 , 2,  3,  or  4)  to  acquire  the  }th,  (j  +  4)^,  (j  +  8)**,  (j  + 
12)fft,  etc.  k-space  lines.  Each  interleaf  is  initiated  with  a  forward 
echo,  b:  A  l-r-l-r-r-l-r-l-etc.  pattern  is  formed  by  acquiring  the  same 
lines  with  each  interleaf  as  in  a,  but  initiating  interleaves  #2  and  #4 
each  with  a  reverse  echo.  Locations  and  relative  intensities  of  ghost 
artifacts  from  these  patterns  are  computed  using  Eq.  [8]  with  h=  4, 
and  are  listed  in  Table  2. 


be  written 


M(x,  y)  =  M{x,  y)  cos  (6{x,  $) 

A/7  2-1  [4] 

+  2  / )  Gir (y  -  /)  sin  (0(x,  /)) 

y=~N/2 


where,  similar  to  that  in  ref.  13,  Gjr(y  -  /)  is  defined  as  the 
Ghost  Kernel  given  by 


N/2-1 


GIr(y  -  /)  s  2  gir(n)  exp  [-i—  n  (y  -  /) J  [5] 


-JV/2 


2tt 


Based  solely  on  the  periodicity  of  glr  with  n,  the  Ghost 
Kernel  can  also  be  written 


GJr(y  -  y)  =  ^  gj^iy-  y  -  [6] 

where  A j=  N/2h  j,j  =  -h+ 1 . h- 1,  and  8  a  delta-like 


function  defined  by 


%'  /  2-it  \ 

8(y)  =  (1/A/)  2  exP  Tr  H 

n=-N/l  \  N  I 

1  ify=  0  Modulo  N 
~  |0  y  integer  ^  0 


[7] 


Eq.  [6]  explicitly  identifies  each  ghost  location  Aj  and 
complex  amplitude  factor  gj.  Substitution  of  Eq.  [7]  intoEq. 
[6],  and  comparison  with  Eq.  [5],  reveals  the  Fourier 
Transform  relationship  between  g}r{n)  and  gy 

gM  =  -i  Eg,  exp  (■* 

[8] 

t,1  I  2ir  \ 

gj  =  0/2 h)  2r  g;,(n)  exp  -j— nAJ 


These  equations  show  that  the  maximum  number  of  ghosts 
is  equal  to  2 h,  and  the  separation  between  the  parent  image 
and  the  two  nearest  (most-overlapping)  ghosts  are  A  {  =  N 
12 h  and  A_j  =  -N/2h  or  ±  1/2 h  the  FOV. 


MATERIALS  AND  METHODS 

Interleaved  gradient-recalled  EPI  was  implemented  on  a 
Signa  Advantage  1.5  T  MR  System  (GE  Medical  Systems, 


FIG.  4.  Image  of  phantom  using  even-number  interleaf  EPI.  In  this 
implementation  of  128  x  128  interleaved  EPI,  each  of  two  interleaves 
acquired  64  k-space  lines,  and  no  lines  were  discarded.  A  /-/-r-r-/-/-r-r- 
etc.  pattern  of  echos  was  formed  by  using  interleaf  #j  (j  =  1  or  2)  to 
acquire  the  Ith,  (j  +  2)m,  (j  +  4)frt,  etc.  k-space  lines.  Each  interleaf 
was  initiated  with  a  forward  echo.  Locations  and  relative  intensities  of 
ghost  artifacts  in  the  image  can  be  computed  using  Eq.  [8]  with  h  =  2. 
The  locations  Aj  and  complex  amplitude  factors  pj  for  these  ghosts, 
denoted  (Aj/A/,  pj),  (with  A/  —  128  image  matrix  size)  are  (-  0), 

(-!4,  Vz\  -  VS,  (Vi,  Vzi  +  Vz).  These  parameters  are  used  in  Eq.  [4]  to 
describe  the  image. 
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FIG.  5.  Axial  images  of  phantom  obtained  from 
single-shot  and  odd-number  interleaf  EPI,  with 
and  without  ghost  artifact  reduction  using  image 
phase  correction,  a:  64  x  64  single-shot  EPI 
with  M-r-l-r-l-r-etc.  pattern  of  echos,  b:  Image 
phase  correction  reduces  average  ghost  inten¬ 
sity  to  2.3%  of  average  parent  image  intensity, 
c:  128  x  128  three-interleaf  EPI,  implementa¬ 
tion  described  in  text  to  generate  M-r-l-r-l-r-etc . 
pattern  of  echos,  d:  Image  phase  correction 
reduces  average  ghost  intensity  to  4.3%  of 
average  parent  image  intensity. 


c 


d 


Waukasha,  WI)  with  Version  5.4  Genesis  software  platform. 
An  insertable  three-axis  gradient  and  radiofrequency  (RF) 
coil  set  (Medical  Advances,  Milwaukee,  WI)  was  used, 
providing  2.0  g/cm  peak  gradient  strength  and  56  psec  rise 
time,  and  permitting  64x64  single-shot  acquisitions  in 
under  40  msec  per  image,  and  with  duty  cycle  permitting 
up  to  9  images/sec.  Figure  1  shows  the  pulse  sequence 
diagram  for  three  interleaf  EPI. 

Specific  interleaved  EPI  scans  were  performed  on  phan¬ 
toms  and  normal  subjects.  Informed  consent  was  obtained 
from  all  subjects.  The  following  parameters  were  common 
to  all  scans:  TE  40  msec,  flip  angle  90°,  FOV  220  mm,  slice 
thickness  6  mm,  slice  gap  2  mm,  and  number  of  locations 
16.  Phantom  studies  were  run  with  TR  3.2  sec  and  subject 
studies  with  TR  2  sec.  Table  1  shows  the  different  inter¬ 
leaved  scans  that  were  performed  and  the  resulting  image 
matrix  sizes.  All  images  were  linearly  interpolated  to 
256  X  256  for  display.  After  image  phase  correction, 
residual  ghost  intensity  was  measured  as  the  average 
intensity  of  the  image  in  the  entire  ghost-only  region  of  the 
image,  divided  by  the  average  intensity  in  the  entire 
parent-only  region  of  the  image,  and  reported  as  a  percent¬ 
age. 

RESULTS 

Predictions  of  Theory 

Typical  interleaved  EPI  data  sets  assembled  with  forward 
and  reverse  echos  grouped  together  produce  ghosts  that 
significantly  overlap  the  parent  image.  Because  image 


phase  correction  requires  at  least  a  one-pixel  wide  strip 
within  the  parent  image  that  does  not  contain  ghost,  the 
technique  cannot  be  applied  with  any  such  interleaved 
data  sets  (unless  the  FOVs  were  made  impracticably  large 
to  separate  the  adjacent  ghosts) .  However,  if  the  interleaved 
EPI  data  could  be  arranged  to  have  purely  alternating 
forward  and  reverse  echos  (i.e.,  1-r-l-r-l-r-l-r-etc.  pattern), 
then  the  resulting  uncorrected  reconstruction  would  have 
only  one  ghost  shifted  by  xh  the  FOV,  and  image  phase 
correction  could  be  applied. 

Odd-Number  Interleaf  EPI 

Even-number  interleaf  EPI  cannot  produce  the  required 
alternating  pattern,  but  odd-number  interleaf  EPI  can.  The 
three  interleaf  EPI  echo  pattern  is  shown  in  Fig.  2.  For 
comparison,  two  different  even-number  interleaf  EPI  echo 
patterns  are  shown  in  Fig.  3.  Figure  3a  shows  a  1-1-7-I-r-r-r-r- 
etc.  pattern  resulting  from  identical  traversals  of  k-space  of 
four  interleaves.  Figure  3b  shows  a  f-r-7-r-r-I-r-i-ete.  pattern 
resulting  from  alternating  the  direction  of  the  traversal  of 
interleaves  #2  and  #4.  More  detailed  consideration  reveals 
that  the  simple  /-r-I-r-I-r-/-r-efc.  pattern  cannot  be  obtained 
from  an  even  number  of  interleaves. 

Experimental  Verification 

Figure  4  shows  a  128  X  128  image  obtained  using  two 
interleaf  EPI  (composed  of  two  128  X  64  interleaves). 
There  are  two  ghosts  shifted  by  ±l/4  of  the  FOV,  arising 
from  the  I-I-r-r-M-r-r-  etc.  pattern  of  phase  errors.  The  ghost 
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FIG.  6.  Sagittal  images  of  head  obtained  from 
single-shot  and  odd-number  interleaf  EPI,  with 
and  without  ghost  artifact  reduction  using  image 
phase  correction,  a:  64  x  64  single-shot  EPI 
with  /-r-/-r-/-r-/-r-efc.  pattern  of  echos,  b:  Image 
phase  correction  reduces  average  ghost  inten¬ 
sity  to  1 .6%  of  average  parent  image  intensity,  c: 
128  x  128  three-interleaf  EPI,  implementation 
described  in  text  to  generate  M-r-l-r-l-r-etc. 
pattern  of  echos,  d:  Image  phase  correction 
reduces  average  ghost  intensity  to  5.6%  of 
average  parent  image  intensity. 


c 


d 


shifted  by  y2  the  FOV  has  zero  amplitude.  Ghost  artifact 
reduction  by  image  phase  correction  is  not  possible,  be¬ 
cause  the  parent  image  is  almost  everywhere  overlapping 
with  ghost.  Figure  5a  shows  an  image  from  64  X  64 
single-shot  EPI,  and  Fig.  5b  shows  removal  of  the  N/2  ghost 
using  image  phase  correction.  Figure  5c  shows  an  image 
from  128  X  128  odd-number  interleaf  EPI  composed  of 
three  128  X  44  interleaves  and  confirms  that  the  sequence 
causes  only  a  single  ghost  shifted  by  l/z  the  FOV.  Figure  5d 
confirms  that  image  phase  correction  largely  suppresses 
this  ghost.  Figure  6  displays  similar  results  obtained  in  the 
sagittal  plane  on  a  normal  subject.  Figure  6a  shows  an 
image  from  64  X  64  single-shot  EPI  with  a  single  ghost 
artifact,  and  Fig.  6b  shows  removal  of  the  ghost  using  image 


phase  correction.  Figure  6c  shows  an  image  from  1 28  X  128 
odd-number  interleaf  EPI  composed  of  three  128  x  44 
interleaves,  and  confirms  that  the  sequence  causes  only  a 
single  ghost  shifted  by  Vz  the  FOV.  Figure  6d  confirms  that 
image  phase  correction  largely  suppresses  this  ghost.  Fig¬ 
ure  7  shows  images  from  256  X  256  odd-number  interleaf 
EPI  (composed  of  nine  256  x  30  interleaves).  These  images 
confirm  that  large  odd-number  interleaf  EPI  produces  only 
a  single  ghost  shifted  by  Vz  the  FOV. 

Comparison  of  Even-  and  Odd-Number  Interleaf  EPI 

Table  3  lists  possible  desired  image  resolutions  and  com¬ 
pares  data  acquisition  time,  number  of  slices,  and  temporal 


FIG.  7.  Axial  phantom  (left)  and  sagittal  head 
(right)  images  obtained  from  odd-number  inter¬ 
leaf  EPI.  256  x  256  combined  data  set  with 
l-r-l-r-l-r-l-r-etc.  pattern  of  echos  was  acquired 
using  nine  256  x  30  interleaves.  Interleaf  #j,  j  = 
1,2, ...,  9,  was  used  to  acquire  the  jth,  (j  +  9),h, 
(j  +  1 8)th,  (j  +  27)^  k-space  line  (number  from  1 
starting  at  top),  and  interleaves  #3,  #5,  #7,  and 
#9  were  initiated  with  reverse  echos  (all  others 
were  initiated  with  forward  echos).  The  first 
acquired  k-space  lines  were  discarded  in  inter¬ 
leaves  #2  through  #9,  as  were  the  last  acquired 
lines  in  interleaves  #1 ,  #5  through  #9.  The  lack 
of  ghosts,  other  than  that  shifted  by  y2  the  FOV, 
validates  the  assumption  that  phase  errors  rep¬ 
resented  by  ±0(x)  are  the  same  in  separately 
acquired  interleaves. 
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Table  2 

Ghost  Complex  Amplitude  Factors  in  h  =  4  Interleaved  EPI  Shown 
in  Fig.  3* 


Ghost 

location 

offset 

f-(-f-/-r-r-r-r-etc. 

pattern 

Ar-Ar-r-/-r-Ae  tc. 
pattern 

-y2 

0 

0 

~3/8 

0.25/-  0.10 

0.25  /+  0.60 

-Vi 

0 

0 

Vs 

0.25  /  -  0.60 

0.25  /+  0.10 

0 

0 

0 

Va 

0.25/ +  0.60 

0.25/-  0.10 

% 

0 

0 

% 

0.25  /  +  0.10 

0.25  /  -  0.60 

*Ghost  position  offsets  are  Ay,  j=-h,  -h  +  7, 1,  as  defined 
in  Eq.  [6],  and  are  listed  here  as  fractions  of  the  FOV  and  relative  to 
the  location  of  parent  image.  Ghost  complex  amplitude  factors  are  gj, 
as  defined  in  Eqs.  [6]  and  [8].  The  echo  pattern  repeats  across  the 
phase  encode  direction,  and  the  amplitude  factors  are  calculated 
assuming  that  the  first  letter  of  the  pattern  gives  the  direction  at  ky  = 
0.  Note  that  both  patterns  generate  ghosts  close  to  the  parent  image 
at  ±1/2/7  =  ±  V8  the  FOV. 

resolution  that  can  be  obtained  using  odd  number  versus 
even-number  interleaf  EPI.  The  table  entries  are  computed 
with  fMRI  neural  activation  experiments  in  mind,  in  which 
temporal  resolution  of  not  worse  than  3000  msec  would 


typically  be  required.  For  a  given  spatial  resolution  and 
number  of  interleaves,  the  right  two  columns  give  the 
number  of  distinct  slice  locations  that  can  be  obtained 
(constrained  by  3000  msec  temporal  resolution),  and  the 
required  TR  (defined  as  the  time  interval  between  acquisi¬ 
tion  of  sequential  interleaves  at  a  given  slice).  The  table 
reveals  that  odd-number  interleaf  EPI  is  complementary  to 
even-number  interleaf  EPI,  and  one  may  argue  that  either 
set  alone  would  be  sufficient  for  the  vast  majority  of  fMRI 
experiments. 

DISCUSSION 

Based  on  our  experience,  reference  scans  are  not  as  reliable 
as  image  phase  correction  (12).  Not  infrequently,  poor 
reference  scans  make  the  image  quality  worse  by  virtue  of 
poor  estimation  of  the  phase  errors  in  the  image  data. 
Reference  scans  will  be  proportionally  more  unreliable  in 
interleaved  EPI  when  separate  phase  corrections  are  re¬ 
quired  for  each  interleaf.  A  recent  comparison  (14)  be¬ 
tween  image  reference  correction  with  two  established 
reference  scan  methods  (5,11)  showed  that  image  phase 
correction  reduced  artifact  more. 

If  the  user  is  willing  to  increase  the  acquisition  time, 
then  the  problem  of  overlapping  ghost  artifact  can  be 
avoided  by  either  1)  doubling  the  FOV  while  maintaining 


Table  3 

Comparison  of  Even-  and  Odd-Number  Interleaf  EPI 


Data  matrix 
(freq  x  ph) 

No.  of 
leafs 

k-Iines/ 

leaf 

Time/leaf 

(ms) 

k-lines 

acquired 

k-lines 

discarded 

No.  of  slices 
with  3000  ms 
temporal  resolution 

Required  TR  (ms) 
for  3000  ms 
temporal  resolution 

64  x  64 

1 

64 

39.9 

64 

0 

27 

3000 

128  X  128 

2 

64 

72.7 

128 

0 

7 

1500 

128  X  128 

3 

44 

50.0 

132 

4 

7 

1000 

128  X  128 

4 

43 

36.4 

128 

0 

7 

750 

128  X  128 

5 

27 

30.7 

135 

7 

7 

600 

192  X  128 

3 

44 

72.5 

132 

4 

4 

1000 

192  x  128 

4 

32 

52.7 

128 

0 

5 

750 

192  x  128 

5 

27 

44.5 

135 

7 

4 

600 

192  x  128 

6 

23 

37.9 

138 

10 

4 

500 

192  X  192 

4 

48 

79.1 

192 

0 

3 

750 

192  X  192 

5 

40 

65.9 

200 

8 

3 

600 

192  X  192 

6 

32 

52.7 

192 

0 

3 

500 

192  X 192 

7 

29 

47.8 

203 

11 

3 

428 

256  X  192 

6 

32 

69.1 

192 

0 

2 

500 

256  X  192 

7 

29 

62.6 

203 

11 

2 

428 

256  X  192 

8 

24 

51.8 

192 

0 

2 

375 

256  X  192 

9 

23 

49.7 

207 

15 

2 

333 

256  X  128 

4 

32 

69.1 

128 

0 

3 

750 

256  x  128 

5 

27 

58.3 

135 

7 

3 

600 

,  256  X  128 

6 

23 

49.7 

138 

10 

3 

500 

256  X  128 

7 

20 

43.2 

140 

12 

3 

428 

256  X  256 

8 

32 

69.1 

256 

0 

1 

375 

256  X  256 

9 

30 

64.8 

270 

14 

1 

333 

256  X  256 

10 

27 

58.3 

270 

14 

1 

300 

256  X  256 

11 

25 

54.0 

275 

19 

1 

272 

Table  lists  only  those  acquisitions  producing  30  ms  <  time/leaf  <  80  ms.  Headings:  Time/leaf:  computed  using  8  us  sampling  rate  and  56  us 
gradient  rise  time.  Data  matrix:  combined  raw  data  array  size,  freq:  frequency  encode  direction,  ph:  phase  encode  direction;  k-lines  discarded: 
Excess  k-space  lines  not  needed  for  reconstruction;  No.  of  slices  . . . :  Estimated  from  time/leaf  and  duty  cycle  limit:  of  local  gradient  coil 
(LGC);  Required  TR  (ms) . . . :  TR  needed  to  acquire  combined  data  set. 
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the  pixel  size  by  doubling  the  number  of  phase  encode 
steps;  2)  retaining  the  FOV  and  acquiring  lines  only  in  one 
direction  with  an  asymmetric  gradient  waveform  as  in 
ABEST  (15,16);  or  3)  doubling  the  number  of  phase  encode 
steps  using  a  standard  readout  waveform  and  performing 
two  reconstructions,  first  using  only  the  forward  lines  and 
then  only  the  reverse  lines,  with  image  averaging  to 
improve  the  signal-to-noise  ratio  (e.g.,  ref.  17  ).  The  odd- 
number  interleaf  EPI  is  valuable  because  in  most  EPI 
applications  it  is  an  optimal  strategy  to  match  the  FOV  to 
the  object  size  and  minimize  the  number  of  phase  encode 
steps,  to  make  image  acquisition  as  rapid  and  with  as  few 
interleaves  as  possible,  while  avoiding  image  wraparound. 
This  strategy  reduces  artifacts  from  physiological  motion 
and  maximizes  temporal  resolution  for  imaging  dynamical 
processes. 
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